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ABSTRACT 

We present a theoretical study of void spins and their correlation properties. 
The concept of the spin angular momentum for an unbound void is introduced 
to quantify the effect of the tidal field on the distribution of matter that make 
up the void. Both the analytical and numerical approaches are used for our 
study. Analytically, we adopt the linear tidal torque model to evaluate the void 
spin-spin and spin-density correlations, assuming that a void forms in the initial 
region where the inertia momentum and the tidal shear tensors are maximally 
uncorrelated with each other. Numerically, we use the Millennium run galaxy 
catalog to find voids and calculate their spin statistics. The numerical results 
turn out to be in excellent agreement with the analytic predictions, both of which 
consistently show that there are strong spatial alignments between the spin axes 
of neighbor voids and strong anti-alignments between the void spin axes and the 
directions to the nearest voids. We expect that our work will provide a deeper 
insight into the origin and properties of voids and the large scale structure. 

Subject headings: cosmology: theory — large-scale structure of universe 



1. INTRODUCTION 

Recent galaxy redshift surveys (e.g., 2dFGRs, Colless et al. 2001) have allowed us to 
map the universe on the largest scale ever and study systematically the detailed properties 
of the cosmic large-scale structures. It is confirmed by these surveys that the galaxies in the 
universe are not distributed evenly but develop pattern of network connected by filaments 
and sheets which in turn are separated by immense voids. 

Voids are defined observationally as large regions with very low galaxy number density 
and shown to occupy approximately 40% of the cosmic volume (Hoyle & Vogeley 2004). The 
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striking vastness of the observed void volumes inspired a flurry of research on the origin and 

properties of cosmic voids (El-Ad & Piran 1997; Schmidt et al. 2001; Peebles 2001; Mathis 
& White 2003; Benson ct al. 2003; Sheth & van de Weygaert 2004; Hoyle & Vogeley 2004; 
Rojas et al. 2004; Colbcrg ct al. 2005; Hoyle et al. 2005; Gottlober et al. 2005; Shandarin et 
al. 2006). It is now generally believed that voids arc originated from the local minima of the 
primordial density fluctuations and expand faster than the rest of the universe. 

The very fact that voids are extremely undcrdense and expand faster has led many 
authors to assume somewhat naively that the shapes of voids should be spherical. In fact 
previous analytical and numerical works on voids arc based on this assumption (e.g., Dubinski 
et al. 1993; Van de Weygaert & Van Kampen 1993; Rriedmann & Piran 2001; Sheth & van 
de Weygaert 2004). 

Very recently, however, Shandarin et al. (2006) investigated systematically the shapes 
of voids by applying the excursion set approach to the density fleld in N-body simulations 
and demonstrated that the intrinsic shapes of voids are far from spherical symmetry. They 
claimed that voids are most vulnerable to external shears due to their low-density, so that 
their shapes are disturbed to be nonspherical in spite of their faster expansion. 

Although Shandarin et al. (2006) focused on the shapes and morphology of voids, we 
note here that if the tidal shear plays a signiflcant role in the case of voids, then it should 
cause not only the shapes of voids to be nonspherical but also voids to acquire the spin 
angular momentum. Our goal here is to construct an analytic model for the spin angular 
momentum of voids originated from the tidal shear effect and test the analytic predictions 
against the data from the recent Millennium run simulation. 

The organization of this paper is as follows. In §2, the analytic model for the void 
spin angular momentum and the void spin correlation properties are presented. In §3, 
the Millennium run galaxy catalog is analyzed to determine void spin statistics, and the 
numerical results are compared with the analytic predictions. In §4, the achievements and 
caveats of our works are discussed, and a flnal conclusion is drawn. 



Let us consider a Lagrangian region in the linear smoothed density fleld. According to 
the linear tidal torque theory (Doroshkevich 1970; White 1984), the spin angular momentum 
of matter which make up this region, J = (Jj), is proportional to the tensor product as 
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THE ANALYTIC MODEL 



Ji ^ ^ijkTjlllki 




-3- 



where I = (lij) is the inertia tensor of the region, and T = (Tij) is the local tidal tensor 
at the center of the mass of the region defined as the second derivative of the perturbation 
potential : = didj^. 

Strictly speaking, equation (1) gives a valid first order approximation to the spin angular 
momentum only when the Lagrangian regions correspond to bound halos. However, to 
quantify the effect of tidal field on voids, we extend the validity of the linear tidal torque 
theory to the case of underdense proto-void regions and define the spin angular momentum 
of a void by equation (1). 

Since a void is an unbound system, its spin angular momentum has a distinct meaning 
from that of a bound halo although defined in a similar way. Basically, the spin angular 
momentum of a void is introduced to quantify how a proto-void region reacts to the effect of 
the local tidal force. The local tidal field will generate non-radial motion of the material in 
the proto-void region, leading the distribution of matter in the void to deviate from spherical 
symmetry. The magnitude of the void spin angular momentum quantifies the degree of this 
deviation, and its direction indicates the direction along which the maximum deviation 
occurs. 

It was shown by N-body simulations that in the overdense proto-halo regions the prin- 
cipal axes of I are in fact strongly correlated with that of T (Lee & Pen 2000; Porciani et al. 
2002). Very recently, Lee (2006) has shown analytically that the degree of the misahgnment 
between the principal axes of I and T decreases monotonically with the linear density of the 
region. In other words, the more underdense a region is, the more misaligned the principal 
axes of I and T are. 

Given these numerical and analytical findings, we propose a new hypothesis that a void 
forms from the initial underdense region where the principal axes of I and T are maximally 
misaligned with each other. This hypothesis combined with equation (1) leads us to expect 
that the alignment between the void spin axes and the principal axes of the local tidal tensor 
should be strongest for the proto-void regions. 

Lee & Pen (2000) showed that the alignment between the spin axis of a Lagrangian 
region and the intermediate principal axes of the local tidal tensor can be best quantified by 
calculating the expectation value of the unit spin given the tidal tensor, and they derived 
an approximate formula for it. Their original formula is characterized by a correlation 
parameter, c, in the range of [0,1], which was introduced to represent the degree of the 
alignment ^ between the spin axes and the principal axes of the tidal tensors. The case of 



In the original formula, Lee & Pen (2000) used a rescaled correlation parameter, a, which is nothing but 
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maximum alignment corresponds to c = 1 while no alignment case to c = 0. For the case of 
halo spins, they found c ~ 0.3 empirically. In accordance of our hypothesis, c = 1 for void 
spins. 

Now, with the value of c set at unity, the formula of Lee & Pen (2000) is rewritten for 
the case of void spins as 

{JiJj\T) = Q + - ^fikfkj, (2) 

where J = J/|J| and T = T/|T| with Tij = Tjj — Tr(T)5jj/3. Equation (2) indicates that the 
spin axes of voids are spatially correlated due to the spatial correlation of the tidal field. The 
alignments between the spin axes of neighbor void spins can be quantified by the two-point 
spin-spin correlation function which was introduced by Pen et al. (2000) for the case of dark 
halos: 

W(r) = (|j(x).j(x + r)r)--. (3) 

Here, the constant 1/3 represents the value of (| J(x) • J(x + r)p) for the case of no alignment. 

For the case of a power-law power spectrum, the void spin-spin correlation function, 
r]v{r), can be approximated (Pen et al. 2000) as 

Here is the correlation function of the linear density field smoothed on a typical La- 
grangian scale of voids, defined as 

sin hv 

^\k)^W\kRv)d\nk, (5) 

where A^(A;) is the dimensionless power spectrum and W{kRv) is the top-hat spherical filter 
of scale radius, i?y . Since voids and galaxy clusters can be thought of the counterparts which 
are supposed to form at the local minima and maxima of the initial density field, respectively, 
we assume the typical Lagrangian scale of voids are same as that of galaxy clusters, setting 
the value of Ry at 8/i~^Mpc. Since the concordance ACDM cosmology is well represented by 
a power-law power spectrum A^(A;) oc A; on the scale of 8/i~^Mpc, one can regard equation 
(3) as a valid approximation. 

It is worth mentioning here that the misalignment between the principal axes of T and 
I (i.e., the value of the correlation parameter c = 1) is characteristic for those underdense 
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regions in the initial density field that will eventually become voids. Although the majority 
of the initial underdense regions will evolve into voids, some underdense regions could merge 
into overdense region to collapse rather than to become voids (Colberg et al. 2005). 

Lee & Pen (2001) showed both analytically and numerically that for the case of dark 
halos the separation vectors between the neighbor halos, r = (n), are ahgned with the 
major axes of the tidal tensors. To quantify the degree of this alignment, Lee & Pen (2001) 
proposed the following formula: 

{hf,\T) = - b Q% - fikfkj^ , (6) 

where f = r/|r| and 6 is a correlation parameter in the range of [0, 1]. The negative sign in 
front of the parameter b indicates that r is aligned with the major principal axis of T unlike 
J. The value of b was determined empirically to be approximately 0.3 for the case of halos 
(Lee & Pen 2001). 

Similarly, we expect that the separation vectors between the neighbor voids are aligned 
with the minor axes of the tidal tensors. Then, equation (6) should be also true for the case 
of voids. We assume further that b has the same value of 0.3. 

Now that the spin axes of voids are aligned with the intermediate axes of the local 
tidal tensors (eq.[2]) while the separation vectors to the nearest voids are ahgned with the 
minor axes of the local tidal tensors, one can expect that the spin axes of voids should be 
anti- aligned with the separation vectors to the nearest voids. 

The expected spin-direction anti-alignment can be quantified by the following correlation 
function (Lee & Pen 2001): 

a;y(r) = (|J-r|^)-l/3. (7) 

It was found by Lee & Pen (2001) that for the case of small-scale galactic halos equation 
(8) can be approximated as RR' /{R^ + i?'^) where R is the typical Lagrangian size of a 
galactic halo and R' = R-\- r. This approximation was made under the assumption that 
the tidal field does not vary significantly within the Lagrangian distance r. For the cases of 
large-scale voids, however, this assumption may not be valid since within the void separation 
the tidal field is likely to vary significantly. Taking into account the spatial variance of the 
tidal field, we find that for the case of voids the spin-direction correlation function is better 
approximated as 
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RESULTS FROM SIMULATIONS 



3.1. Identifying Voids 



To test the analytic predictions made in §2, we use the Millennium run galaxy catalog 
carried out by the Virgo Consortium (Springel et al. 2005). The catalog consists of total 
8964936 galaxies at redshift z — in a, periodic box of linear size 500/i~^Mpc, with infor- 
mation on various galaxy properties such as position, velocity, mass, magnitude, and star 
formation rate. 

Voids are identified in the catalog by the void-finding algorithm described in Hoyle & 
Vogeley (2002, hereafter HV02) which was in fact first suggested by El- Ad & Piran (1997). 
To apply the HV02 algorithm to the Millennium run catalog, we first determine the values 
of the key parameters, / and s^. The two parameters, / and Sc, represent the distance 
within which a field galaxy is required not to have any three neighbors and a threshold for 
the minimum linear size of a void, respectively. For the Millennium run catalog, the key 
parameter values we obtain are I — 2.44/i~^Mpc and Sc — 6/i~^Mpc, respectively. For the 
detailed description of the procedures to identify voids from the Millennium run catalog, see 
our companion paper (D. Park & J. Lee 2006 in preparation). 

The total 24037 voids are found from the catalog. The volume Ty of each void is 
measured by means of the Monte-Carlo integration method described in HV02. The effective 
radius of each void is then determined as i?e = (3ry/47r)''^/^. The mean value of Re is obtained 
to be Re = 10.45/i-^Mpc. 

The density of each void, dy, is measured by counting the number of void galaxies: 
Sy = {ny — fig) /fig where ny = Ny/Ty is the number density of void galaxies, and fig is the 
mean number density of the total catalog galaxies. Figure 1 plots the probability distribution 
of 5v- As can be seen, the voids found in the catalog are extremely underdense with the 
mean density 6v — —0.92. 



We first choose only those voids which contain more than 30 void galaxies. Total number 
of such voids is found to be 13507. The specific spin angular momentum, j, of each selected 
void is determined as 



3.2. Measuring the Void Spin Angular Momentum 




(9) 
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where and My = X^a^ the mass of the ct-th galaxy and the total mass of all 

galaxies in each void, respectively. The position and the velocity of the a-th void galaxy in 
the center of mass frame are represented by and v^, respectively. We rescale this specific 
spin angular momentum to be dimensionless as 

J = , (10) 

V2MvReV 

with V'^ = GMy/Re-, which is analogous to the definition of the halo spin parameter (Bullock 
et al. 2001). Figure 2 illustrates a few examples of the voids identified from the catalog with 
the directions of the void spin axes marked as short solid line on the voids. 

Figure 3 plots the distribution of j for the voids from the catalog as histograms, which 
shows that p{j) has a log-normal shape. We fit the following log-normal formula to the 
histogram with adjusting the values of crj and Jq: 

The best-fit values of cxj and jo ai'e found to be 0.72 and 0.91, respectively. As can be seen, 
the value of j spreads over a large range, revealing the strong shear effect on voids. 



3.3. Determining the Void Spin-Spin and Spin-Density Correlations 

We measure the cosine of the relative angle, cos 9, between the spin axes of each void pair 
in the catalog, and determine the average of cos^ as a function of the separation distance 
between the two voids. For this, we also use only those voids which have more than 30 
galaxies. The two-point spin-spin correlation function (eq.[3]) is then determined simply by 
subtracting one third from the average, (cos^ 9) . 

Figure 4 plots the numerical result as solid squares with errors. The errors are calculated 
as the standard deviation of cos^ 9 for the case of no correlation. The horizontal dotted hue 
corresponds to the case of no alignment. As can be seen, there is a clear signal of alignment 
up to the distance of 20/i~^Mpc, which is almost an order of magnitude stronger compared 
with the case of halos. 

The analytic prediction (eq. [4]) is also plotted as solid line for comparison. As can be 
seen, the analytic prediction is in excellent agreement with the numerical result. It is worth 
emphasizing that the analytic prediction is not a fitting model but derived from first physical 
principles under the assumption that voids form in the initial regions where the tidal and the 
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inertia tensors are maximally uncorrelated with each other. The good agreement between 
the numerical result and the analytical prediction supports this idea. 

In a similar manner, we also determine the void spin-direction correlation function 
(eq.[7]). Figure 5 plots the numerical result as solid squares with errors. As can be seen, 
there is a clear signal of anti- alignment between the void spin axes and the direction toward 
the nearest neighbor voids, which is also an order of magnitude stronger compared with the 
case of halos. The analytic prediction (eq.[8]) is also plotted as solid line for comparison. We 
find a good agreement between the analytical and the numerical results for the spin-direction 
correlation, too. 

4. DISCUSSION AND CONCLUSION 

In the light of the work of Shandarin et al. (2006) who claimed that the tidal shear must 
play a strong role in case of voids, we have proposed a scenario that voids are originated from 
the regions in the initial density field whose inertia tensors are completely uncorrelated with 
the local tidal tensors. Since the alignment between the principal axes of the inertia and 
tidal tensors plays a role of reducing the tidal shear effect according to the linear tidal torque 
theory, the hypothesis implies that voids should have very high spin angular momentum. We 
define the concept of the spin angular momentum for an unbound void to quantify the effect 
of the tidal field on the distribution of matter that make up voids. 

Our analytic model based on the linear tidal torque theory have predicted that there is a 
strong spatial alignment between the spin axes of neighbor voids and a strong anti- alignment 
between the spin axes of voids and the directions to the nearest voids. 

Our predictions have been tested against the Millennium run galaxy catalog where 13507 
large voids with more than 30 galaxies are identified. We have measured numerically the 
distribution and correlation functions of void spins, and found that the predictions of our 
analytic model are in excellent agreement with the numerical results. 

The success of our work is, however, subject to a couple of caveats. The first caveat 
comes from our implicit assumption that the matter voids are same as the galaxy voids. In 
our analytic model, we consider the matter voids. While in the numerical analysis we deal 
with the galaxy voids. Although the two types of voids are not necessarily same, they should 
be at least closely correlated with each other. 

The second caveat hes in the fact that there is no standard way to identify voids. Unlike 
the case of halo-finding, the use of different algorithms could produce different results for 
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some properties of voids. Especially, the shapes and volumes of voids can differ significantly 
by the underlying criterion of the void-finding algorithm. If one use an algorithm which 
regards a void as a spherical region of low density, then all voids found by that algorithm 
would be spherical. In spite of this ambiguity in finding voids, we expect that our result on 
the correlation properties of the void spin angular momentum would not change significantly 
by the different choice of the void finding algorithm, since the correlation properties are 
basically determined not by the shapes of voids but by the galaxies in voids. 

Another issue we would like to discuss here is a possibility of testing numerically our 
key assumption that a void forms in the initial undcrdcnsc region where the misalignment 
between the principal axes of the inertia and the tidal tensors becomes maximum. Although 
Lee (2006) has already shown analytically that our assumption is true, it would be highly 
desirable to test it against numerical data from the Millennium run simulation. For this, 
it would be necessary to find the Lagrangian positions of all the particles that make up 
each void and to determine the inertia momentum tensor of each Lagrangian proto-void 
region in the principal axis frame of the initial tidal field. As the initial density field of 
the Millennium run simulation is not available to the public yet, the numerical test of our 
assumption is postponed. We hope to report the test result elsewhere in the future. 

Finally, we conclude that our work on void spins will provide a deeper insight into the 
origin and properties of the large scale structure in the universe. 

We thank the anonymous referee for helpful suggestions. The Millennium Run simu- 
lation used in this paper was carried out by the Virgo Supercomputing Consortium at the 
Computing Center of the Max-Planck Society in Garching. The semi-analytic galaxy cat- 
alogue is publicly available at http://www.mpa-garching.mpg.de/galform/agnpaper. This 
work is supported by the research grant No. ROl-2005-000-10610-0 from the Basic Research 
Program of the Korea Science and Engineering Foundation. 
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Fig. 1. — Void density distribution measured from the Millennium run galaxy catalog. 
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Fig. 2. — Examples of the voids identified in the Millennium run galaxy catalog. The straight 
line on each void represents the spin axis. 
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Fig. 3. — Distribution of the rescaled spin angular momentum of voids. The histogram 
represents the result from the Millennium run simulation while the solid line is the log- 
normal fit (eq.[ll]). 
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Fig. 4. — Void spin-spin correlation function (eq.[3]). The solid squares with errors represent 
the result from the Millennium run simulation while the solid line is the analytic prediction 
(eq.[4]). The horizontal dotted line corresponds to the case of no alignment. 
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Fig. 5. — Void spin-density correlation function (eq.[7]). The solid squares with errors 
represent the result from the Millennium run simulation while the solid line is the analytic 
prediction (eq.[8]). The horizontal dotted line corresponds to the case of no ahgnment. 



